Spin-accumulation in small ferromagnetic double barrier junctions 
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The non-equilibrium spin accumulation in ferromagnetic 
double barrier junctions is shown to govern the transport in 
small structures. Transport properties of such systems are 
described by a generalization of the theory of the Coulomb 
blockade. The spin accumulation enhances the magnetore- 
sistance. The transient non-linear transport properties are 
predicted to provide a unique experimental evidence of the 
spin-accumulation in the form of a reversed current on time 
scales of the order of the spin-flip relaxation time. 
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In the seventies it was understood that electron trans- 
port in tunneling and heterostructures involving meta^i 
ferromagnets is associated with non-equilibrium spins 
Compared to other time scales in electron transport the 
spin relaxation time is generally very long at low tem- 
peratures, being limited only by scattering at paramag- 
netic impurities and by spin-orbit scattering. The spin- 
relaxation time and the spin-diffusion length which gov- 
ern the spin accumulation has been measured by John- 
son in polycrystalline gold films H The concept of non- 
equilibrium spin accumulation plays an important role 
in the Boltzmann theory of transport of the giant mag- 
netoresistance in the current perpendicular to the plane 
(CPP) configuration na However, the experimental evi- 
dence for the spin accumulation is indirect at best. It can 
be shown that in the linear response regime the spin- and 
charge- distribution functions can be completely inte- 
grated out of the transport problem, which then depends 
exclusively on the scattering probabilities and the applied 
bias .El In this Rapid Communication we show theoreti- 
cally how unambiguous evidence for a non-equilibrium 
spin accumulation can be obtained by the DC and AC 
response of ferromagnetic double barrier junctions in the 
non-linear regime. These junctions have to be small in 
order to observe big effects, which means that the com- 
plications of the Coulomb blockade have to be taken into 
account (for a review see Ref. [ |(| ) . To this end we have 
to extend very recent theories of the Coulomb blockade in 
ferromagnetic double barrier junctions!] to inclwie time 
dependence and a non-zero spin relaxation timell Ono et 
al. succeeded in fabricating a ferromagnetic single elec- 
tron transistor,EI which in principle can be used to test 



our predictions. Coulomb charging-effects have also been 
seen in discontinuous multilayeral3 and in small cobalt 
clusters. til 

We first show that the spin accumulation in a ferro- 
magnetic double barrier junctions becomes relevant when 
the number of electrons in the island between the tun- 
neling barriers is relatively small. In ferromagnetic struc- 
tures where the tunneling rates depend on the electron 
spin, a finite current through the system is accompanied 
by a spin current out of or into the island (ds/dt) tr . This 
creates a non-equilibrium excess spin s on the island, 
which decays with the spin-flip relaxation time r s f so that 
in steady state (ds/dt) tr = s/t s {. Energy relaxation is 
much faster than spin relaxation, so that the occupation 
of the states for each spin direction can be described by 
Fermi distributions Jj The non-equilibrium spin accumu- 
lation on the island is equivalent to a chemical potential 
difference A/i between the spin up and the spin down 
states. Since spin relaxation is slow and the structures 
of interest are small, A/j, is uniform over a sufficiently 
small island. In terms of the typical single-particle en- 
ergy spacing (or inverse energy density of states at the 
Fermi energy) 5 we have A/i = sS. Spin accumulation 
may be expected to interfere with the transport proper- 
ties when A/x is of the same order as the applied volt- 
age V . The spin current is of the same order as the 
current, e(ds/dt) tI ~ I ~ V/R, where R is the typical 
junction resistance. The non-equilibrium spin accumula- 
tion is therefore important when the spin-relaxation time 
and/or the single-particle energy spacing are sufficiently 
large: 



TssS/h > R/Rk: 



(1) 



where the quantum resistance is Rk = h/e 2 . The spin 
flip relaxation time in polycrystalline aluminium is r s f 
10 10 sOl (lCT 8 s in single-crystal aluminium at T = 4.3-fd 
for gold a The single-particle energ 



) and r s f ~ lCP^s tor goldJJ Trie smgle-particle energy 
spacing on the island is roughly 6 ~ Ep/N, where N is 
the number of atoms on the island and Ep ~ \§eV is the 
Fermi energy. In an Al island with less than 10 6 atoms 
(10 s atoms in single-crystals) the spin accumulation may 
therefore be expected to play a significant role. "Mod- 
ern" metals, like arm-chair nanptubeslij or (magnetic) 
semiconductor heterostructurea!3 can also be interesting 
as island materials. The first because of a possible huge 



1 



spin-flip relaxation time and the latter since islands con- 
taining a small number of electrons can be created by 
depletion of the two-dimensional electron gas .19 

In small systems where Eq. ([!]) is satisfied the spin- 
flip relaxation time is longer than the charge relaxation 
time RC (C is the capacitance of the island). This can 
be seen from Eq. (0), r s f > (2E C /S) ■ RC, and not- 
ing that the charging energy is larger than the single- 
particle energy spacing except in few-electron systems, 
E c /S ~ (e 2 /E F a) N 2 ' z (e 2 /E F a ~ 1). Hence the long- 
time response of the system is dominated by the spin 
dynamics. 

We consider a normal metal island attached to two 
ferromagnetic leads by two tunnel junctions. We assume 
collinear magnetizations in the leads and disregard size 
quantization. The tunnel junctions are characterized by 
a capacitance Cj and magnetic configuration-dependent 
conductances Gi a , where i = 1, 2 denotes the first and the 
second junction and tr denotes up (+) or down (— ) spin 
electrons on the island. There is a source-drain voltage 
V between the right and the left reservoir and a gate 
voltage source coupled capacitively to the island. Here 
we consider the situation with a maximum Coulomb gap 
wherg the offset charge controlled by the gate voltage is 
zero.B 

We proceed from the assumptions of the orthodox 
theory, i.e. Gi G < Gk neglecting co-tunneling,EJ with 
the difference that the transition rates becomes spin- 
dependent. The transition rate from the left reservoir 
to the island is 



\+i,n = -^G la F{E 1 {V,q)-aAfi/2), 



(2) 



where the energy difference associated with the tunnel-, 
ing of one electron into the island through junction i isu 
Ei(V,q) = n.eV + e(q - e/2)/(Ci + C 2 ), the charge on 
the island is q = —ne, the total capacitance is 1/C — 
1/Ct + 1/C 2 , Ki = C/Ci, F{E) = E/[l - cxp(-E/k B T)} 
and fc_eT is the thermal energy. The spin balance is 



ds 
dt 



(3) 
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where the spin- relaxation rate is (ds/di) , = — s/r s f = 
— A/i/c)T s f, r s f is the spin-flip relaxation time and S^ 1 is 
the density of states at the Fermi level in the island. The 
spin balance (Q) can be written in the stationary case as 
I s = e(ds/dt) tr — G s 2A/i/e, where the "spin relaxation 
conductance" is introduced as G s = e 2 /2Sr s f. The master 
equation^ determines the probability p n to have n excess 
electrons on the island. The current through the first 
junction is I 1 = (ij +if ), where the current of electrons 

with spin a is I{ = e^„Pn(r 1<T „+i,„ - r lcr „_i >n ) and 
there is a similar expression for the current through the 
second junction I 2 = {l\ + 1^)- The spin current is 
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(4) 



In the Coulomb blockade regime the current is zero, 
1 = 0, And it can be shown that A^i vanishes, as 
expected.tj The Coulomb gap in the low-temperature 
current-voltage characteristics is thus not modified by 
the non-equilibrium spin accumulation. We also want 
to point out that for symmetric tunneling junctions 
Gii /Gii = G21/G21 the non-equilibrium spin accumu- 
lation vanishes and our theory reduces to those in Refs. 

[§■ 

In this orthodox model the problem can be mapped on 
the equivalent circuit in Fig. |l| by introducing the "spin 
capacitance" C s = e 2 /25, so that 

(es)/2 = C s (A/i/e), A^/s = e 2 /{2C s ) = S. 

This "charging energy" of the spin capacitance is thus 
simply the single-particle energy cost of a spin-flip, 5, or 
more generally, the inverse of the magnetic susceptibility 

We solve the general problem for the steady state as 
well as for the time-dependent properties by numerically 
integrating the master equation and the spin balance, Eq. 
(H). We choose symmetric capacitances Ci = C2 = C 
in our calculations. Thus the important energy scale 
is the Coulomb energy E c = e 2 /2C and the other rel- 
evant energies are renormalized by E c . The thermal 
energy is ksT = 0.05E C . The spin-dependent junction 
conductances are described in units of the average junc- 
tion conductance G and the currents are normalized by 
Ge/2C. In the parallel configuration, the conductances 
are G\ a = d(l + <jP)/2 and G\ a = G 2 (l + aP)/2, 
where P is the polarization of the ferromagnets. In 
the antiparallel configuration G^ = Gi(l + crP)/2 and 
G£ = G 2 (l-aP)/2. 

We consider first the steady state transport proper- 
ties where the spin capacitance C s does not contribute. 
The junction magnetoresistance is the relative difference 
in the resistance when switching from the antiparallel to 
the parallel configuration. In the absence of the non- 
equilibrium spin accumulation, the junction magnetore- 
sistance vanishes for the F/N/F junction. The spin accu- 
mulation causes a non-zero magnetoresistance. We show 
in Fig. ^| the calculated junction magnetoresistance for 
Gi/G = 1, G 2 /G = 2 and a polarization P = 0.4 in 
the limit of slow spin-relaxation G s /G = (upper curve) 
and fast spin-relaxation G s /G = 5 (lower curve). We see 
the magnetoresistance oscillations as a function of the 
source-drain voltage.13 The amplitude of the oscillations 
decreases with increasing source-drain voltage, where the 
Coulomb charging is less important .0 The period of the 
oscillations is close to 2E C for our system. There is only 
a small distortion of the shape of the magnetoresistance 
oscillations with increasing spin-relaxation rate in the is- 
land. The magnetoresistance and its oscillations are no- 
ticeable even when the spin-relaxation conductance is of 
the same order as the tunnel conductances, in agreement 
with Eq. (Q). In the absence of the Coulomb charging 
energy, the tunnel magnetoresistance is 
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TMR = P z 



1 - P 2 7 2 



(5) 



where 7 = (Gi — G 2 )/(Gi + G 2 ) is a measure of the asym- 
metry of the junction conductances and a 2 — 4G s /(Gi + 
G2) determines the reduction of the magnetoresistance 
due to the spin-relaxation. For a high source-drain bias 
when the Coulomb charging effects are negligible, the nu- 
merical results agree well with Eq. (|), TMR= 11% for 
G s /G = and TMR= 2% for G s /G = 5. 

For the transient response in the antiparallel config- 
uration we use P = 0.5, G x /G = 1.3, G 2 /G = 2.6 
and G s /G = 0.3. Let us consider first a fixed source- 
drain voltage at a high bias until the system is stationary 
and then lower the source-drain voltage. We have used 
T sf = imC (e.g. E c = 0.2 meV and R/R K = 10 gives 
RC = 2 • 10~ n s). The initial high bias is V t = IQE C 
which gives a stationary current of 2, = 6.2Ge/2C and 
we investigate the behavior of the transient current when 
the final source-drain bias is below, Vf = (// = 0), 
and above the Coulomb charging energy, Vf = 4E C 
(If = 2.1Ge/2G). We show in the upper panel in Fig. || 
the current through the first and the second junction for 
Vf — 4:E C (upper curves) and Vf = (lower curves) after 
the source-drain voltage is changed at t = 0. It is clearly 
seen that the relaxation of the current is slow on the 
time scale RC. For time scales less than RC, we see that 
the current through the first and the second junction are 
not the same due to the charge depopulation in the is- 
land. The average of the upper curves (Vf = 4E C ) where 
the final source-drain voltage is well above the Coulomb 
blockade energy, follows to within 10-20% the description 
given by the equivalent circuit neglecting the Coulomb 
charging effects described below (|J), (^) and (||) accord- 
ing to which the spin accumulation time is r sp i n = 2ARC. 
When the source-drain voltage is switched off (Vf = 0), 
we see that the transient current is negative. However, 
the spin accumulation time is much longer in this case, 
T S pi n ~ r s f = 10i?G. This discrepancy becomes more 
evident when we consider the relative change of A/i or s: 



D(t) 



s(t = 00) — s(t) 



?(t = 00)- s(t = 0) 



In the lower panel of Fig. y we show the calculated time- 
dependent relative change D(t) in the situations Vf — 
AE C (upper solid curve) and Vf — (lower solid curve), 
which are found to be remarkably different. 

In order to understand the dynamics it is useful to 
inspect the device without the Coulomb charging effects, 
i.e. the capacitances C± and G 2 in the equivalent electric 
circuit in Fig. El We set the voltage on the left lead to zero 
and apply a time dependent potential V(t) to the right 
lead. The complex impedance Z sp i n (u>) = V(uj)/I(uj) is 



G\G 2 G\^G 2 [ — GiiG 2 -\ A/j,(uj) 



Gi + G 2 



(Gi + G 2 ) eV(w) 



(6) 



A/x(cj) _ 1 Gi T G 2 x - GijG 2 | 

eV(w) ~ 1 + twTiph, (G s + G')(Gi+G 2 )' 

Here the spin accumulation time is 



' spin 



G s + G' ' 



(7) 



(8) 



where 



where 1/G' = l/(Gi T + G 2T ) + 1/(G U + G 2i ). From 
the relations (^) and (0) we see why switching-off the 
source-drain voltage (Vf = 0) reverses the transient cur- 
rent as found in the upper panel in Fig|j. Without the 
Coulomb blockade this transient decays on the time scale 
T sp i n . In the limit that the junction conductances are 
much smaller than the spin conductance, the spin ac- 
cumulation time (||) reduces to the spin-flip relaxation 
time, r S pi n ~ T s f . In the opposite limit where the junc- 
tion conductances are much larger than the spin con- 
ductance, the spin accumulation time is r sp i n ~ C S R. 
The spin-capacitance is much larger than the charge- 
capacitance G in the regime where the orthodox theory 
is valid (5 <g; Ec) and thus the spin accumulation time 
is much larger than the charge-relaxation time. 

The dashed lines in the lower panel in Fig. || corre- 
spond to the spin accumulation time in the absence of 
charging, r sp i n = 2.4i?G as well as to the spin-flip relax- 
ation time T s f = 10i?G. We see that the calculated spin 
accumulation time agrees well with the equivalent circuit 
described above (0) for Vf — 4i? c , but disagrees with 
this expression for Vf = where the spin accumulation 
time is close to r s f. The latter is a result of the Coulomb 
charging which is seen to affect the spin accumulation 
time. In this case the non-equilibrium spin accumulation 
decays slower since the spins must relax through the spin- 
conductance G s on the island and the transport through 
the junctions is suppressed. In this situation, the relax- 
ation time of the non-equilibrium spins and the current 
for long times is equal to the spin-flip relaxation time T s f , 
as observed in Fig. [|. 

It should be noted that the magnon assisted inelastic 
tunneling, which reduces the TMR, gives negligibly small 
contribution in our case because of a magnon excitation 
gap, pijesumably due to magnetic anisotropy and/or size 
effects.ta This magnon gap is larger than the bias volt- 
age applied in our study. For verysmall islands like the 
metallic cluster studied in Ref. [ [TlJ] when the Coulomb 
charging energy is larger than the magnon gap, magnon 
inelastic tunneling can interfere with the Coulomb charg- 
ing effects. 

In conclusion, we have investigated the influence of 
a non-equilibrium spin accumulation on the transport 
properties of a ferromagnetic single-electron transistor. 
For a F/N/F junction we find a finite magnetoresistance 
due to the non-equilibrium spin accumulation. The spin 
accumulation can have a drastic effect on the transient 
transport properties. A transient response can be found 
on time scales much larger than the charge relaxation 
time RC. The same slow response is also expected if 
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other external parameters such as the gate voltage, or 
the magnetization are changed. 
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FIG. 2. The junction magnetoresistance in the limit of no 
spin-relaxation in the island (Gs /G = 0) and fast spin relax- 
ation (Gs/G = 5). 

FIG. 3. The current as a function of time (upper panel). 
The relative change of the non-equilibrium spin as a function 
of time (lower panel). The source-drain voltage is switched 
from Vi = 1QE C to V f = 4E C or V f = at t = 0. 
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FIG. 1. The equivalent circuit for the current-voltage re- 
sponse of the system. 
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